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ABSTRACT 

Aims. The standard thin accretion disk model has been successfully used to explain the soft X-ray spectra of Galactic black hole 
systems and perhaps the UV emission of Active Galactic Nuclei. However, radiation pressure dominated disks are known to be 
viscously unstable and should produce large amplitude oscillations that are typically not observed. Instead, these sources exhibit 
stochastic variability which may naturally arise due to viscous fluctuations in a turbulent disk. Here we investigate whether these 
aperiodic viscous fluctuations can stabilize the inner radiation pressure dominated disks and hence maybe the answer to a forty year 
old problem in accretion disk theory. 

Methods. The structure and evolution of a time dependent accretion disk around a black hole is solved numerically. We incorporate 
fluctuations in the disk by considering stochastic variations in the viscous parameter a on the local viscous time-scale. We study both 
locally stable disks where the viscous stress scales with the gas pressure and the locally unstable disks where the stress scales as the 
total pressure. We consider steady state disk parameters pertaining to both stellar and super-massive black holes. 
Results. For locally stable disks, the power spectra of the luminosity variations are found to inversely scale with frequency i.e. 
P{f) oc 1// as expected from analytical studies of linear viscous fluctuations. For unstable disks, where the viscous stress scales with 
the total pressure, the standard oscillatory solutions are seen when the viscous fluctuation amplitude is small. Increasing the fluctuation 
amplitude decreases the amplitude of these oscillations till for a sufficiently large (near unity) fluctuation amplitude, the oscillatory 
behavior disappears and the luminosity variation of the disk becomes stochastic. This study may explain why many accreting black 
holes, despite having accretion rates large enough to develop the radiation pressure instability, do not exhibit large amplitude regular 
outbursts, or like the micro-quasar GRS 1915-H105, produce diminished versions of them only in some particular spectral states. 
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(^ • 1 . Introduction putation based on the standard theory, but incorporating general 

T— I ' relativistic effects and detailed radiative transfer have been used 

(3 ; The standard accretion disk model (Shakura & Sunyaev 1973) successfully used to fit the soft component. Such analysis can 

. , has been the basic framework by which the accretion process constrain the mass and spin of the black holes and opens up the 

^ around black holes is understood. The model, which is also re- jbility of testing General Relativity in the strong field limit. 

^ ferred to as the a-disk model, assumes a thin Keplerian accre- ^j^^^^ ^j^^j^ ^j^^^.^ ^^ ^^^^ ambiguity regarding the origin of 

^ tion disk where the unknown turbulent viscosity is parametrized ^j^^ j^^^^ ^-ray emission from black hole binaries, there is near 

^ such that the VISCOUS torque is proportional to the total pres- consensus that the soft component is due to a standard accretion 

• ^ !r/f ■ ^'' : "'^'r ^^ niodel predicts that for super-massive ^-^^^^ ^^-^^ ^jj^^^ ^^^ ^ ^^^^^^^ ^^^ unprecedented understand- 

k^ black holes the peak radiation will be emitted in the ultra-violet • ^ these sources 

'^^ and this is often identified with the "Big Blue Bump" in Active 
C^ : Galactic Nuclei (AGN). For stellar mass black holes, the emis- There has been a long standing concern regarding the stan- 

- - - sion peaks in soft X-rays and such a thermal soft emission is dard accretion disk model, that the disk is both thermally and vis- 

indeed observed in Galactic black hole binaries. cously unstable when radiation pressure is dominant (Lightman 

Galactic black hole binaries exhibit different spectral states & Eardley 1974). This occurs in the inner regions of the disk 

(see e.g., for a review. Done et al. 2007). In the hard/low state the ^0^"^ ^ %^ solar mass black hole, when its luminosity exceeds 

spectra is dominated by an hard X-ray power-law that extends to ^3x10^^ ergs/s. X-ray binaries can typically exceed this lu- 

~ 100 keV and a weaker soft black body like component. The mmosity which is substantially smaller than than the Eddington 

soft component is generally modeled as a truncated standard ac- 1™!^ of ~ 10^^ ergs/s. For AGNs, harboring lO^M^ black holes 

cretion disk and the primary hard power-law is believed to arise *e disk should be unstable for luminosities exceeding > 5 x 10^^ 

from a hot inner disk (e.g., Shapiro, Lightman & Eardley 1976). ^rgs/s sigmficandy less than the Eddington limit of ~ 10^^ ergs/s. 

In the soft state, the spectrum is dominated by a thermal emis- Thus, for typical black hole binaries and AGNs, the standard ac- 

sion from the standard accretion disk. For X-ray novae, lumi- cretion disk is unstable and not hence not applicable, 
nosity changes by several orders of magnitude, yet the inferred The extension of the radiation pressure dominated part of 

inner radius of the standard accretion disk remains constant, as the disk depends on the accretion rate, the black hole mass and 

expected from the standard disk theory. In fact, spectral com- the viscosity. The parameter regions of the radiation pressure 
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dominated disks in Galactic black holes and AGN were shown 
recently in Janiuk & Czemy (201 1); the approximate fitting for- 
mula for this dependence is: 



log — ^ = log m + 2.5 +0.1 log 

^ Rs 0.01 



(1) 



for a Galactic black hole of 10 M© at accretion rates above 
log merit = -1.5 where m - M/Mndd and Rs - 2GM/c^. For 
a supermassive black hole of lO^M© at accretion rates above 
log merit - -2.5, the radial range is: 



log ^ = 0.5 logm + 2.7 + 0.1 log ^ 



(2) 



Initially it was postulated that the unstable inner regions of 
the standard disk may get converted into a hot two temperature 
region (Shapiro, Lightman & Eardley 1976). The hot inner disk 
could then be the source of the hard X-ray emission and is a con- 
sistent model for Galactic black holes in the hard state. However, 
in the soft state, where the luminosity could be even higher than 
the hard one, there is a soft thermal-like emission indicating that 
the standard disk somehow survives the instability caused by the 
radiation pressure instability. Moreover, time-dependent numer- 
ical evolution of the unstable disk, shows that instead of the disk 
getting transformed to a hot one, it undergoes cyclic behavior, 
which result in large amplitude oscillations in luminosity. These 
large variations in luminosity have not been observed. 

Attempts have been made to stabilize the disk, or at least 
reduce the destabilizing effect, by accounting for the fraction 
of power that maybe transferred to an overlying corona or jet. 
This indeed increases the critical accretion rate, above which the 
instability occurs by a factor of few, nevertheless the instabil- 
ity persists and the model predicts large amplitude oscillations 
at high luminosities. Further, assuming that there is a vertical 
height dependence on the energy dissipated in the disk, where 
a significant fraction of the energy is dissipated near the pho- 
tosphere and not the mid-plane, also does not completely sta- 
bilize the disk. Finally, one can bring to question, the form of 
the viscous stress chosen in the standard model, where the stress 
is taken to be proportional to the total pressure, Tr^ = aPiot- 
Indeed, if the viscous stress scales with only the gas pressure 
Tr^ = af gas, the disk would be unconditionally stable. However, 
it is not clear why such an ad hoc prescription would be a 
physical one. If the turbulence length scales with the height of 
the disk, but the eddy speed scales with the gas sound speed, 
Cs = y/yPgds/p (instead of the total sound speed c, - ^J y P tot Ip), 
the viscous stress would be oc -y/Pgas^tot- However, even for such 
a prescription the disk is not stable, but undergoes luminosity os- 
cillations with diminished amplitude. A combination of such a 
viscous prescription and assuming that a fraction of the power 
is transferred to a corona, may explain the variability of the 
micro-quasar GRS 1915 + 105, in certain spectral states (Janiuk 
& Czerny 2011). GRS 1915+105 is not typical and the appar- 
ent stability of other bright X-ray binaries and even for GRS 
1915+105 in its stable states, remains a puzzle (see e.g, Neilsen 
etal.2011). 

Instead of showing large amplitude oscillatory behavior, 
black hole binaries and AGNs exhibit stochastic variability with 
a high frequency cutoff (Nowak et al. 1999). In the hard state 
of the X-ray sources, the photons from the disk are being up- 
scattered due to the Comptonization in the hot plasma. The dis- 
persion of the scattering times may lead to the attenuation of 
the observed variability amplitude. Also, the dependence of the 
power density spectrum (PDS) with energy may be due to direct 



fluctuations in the Comptonizing medium (Nowak et al. 1999). 
In the soft state, the observed radiation is mainly from the ac- 
cretion disk, whose dynamical variability may be coupled to the 
emission processes in a non-linear way. Also, the effect of re- 
flection of the hard X-rays from the disk surface may lead to a 
feedback effect and influence the variability amplitude (Uttley & 
McHardy 2005). For some sources, gravitational light bending 
may also be affecting the variability of the source (NGC-6-30- 
15; see Miniutti & Fabian 2004). In any case, the amplitude of 
the observed variability in both Galactic black holes and AGN is 
of the order of 30%. 

One of the models for the physical origin of the disk vari- 
ability is stochastically driven damped harmonic oscillations at 
the local epicyclic frequency (Misra & Zdziarski 2008). 

Alternatively or additionally, the variability may be due to 
variations in the local viscous stress, which induces local accre- 
tion rate fluctuations, which in turn induces and adds to variabil- 
ity downstream at smaller radii (Lyubarskii 1997). The model 
predicts a 1 // like variability for the inner disk accretion rate and 
a cutoff at the frequency corresponding to the viscous time-scale 
of the inner region. Indeed, the observed power spectra of black 
hole binaries and AGNs can be explained in terms of this viscous 
fluctuation models (e.g., Churazov et al.2001; Uttley et al.2005; 
Titarchuk et. al. 2007, Ingram & Done 201 1). More interestingly, 
since in this model, the disturbance travels from the outer to the 
inner regions and if higher energy photons are preferential ly pro- 
duced in the inner regions, there would be time-lag between the 
high and low energy photons. The cylindrical geometry would 
make this time-lag inversely proportional to frequency (Misra 
2000). The observed time-lags in black hole binaries and AGNs 
has this frequency dependence (Nowak et al. 1999) which gives 
strong credibility to the viscous fluctuation model, which other- 
wise would be difficult to explain. Theoretically too, one expects 
that a turbulent viscous stress, perhaps arising from the magneto- 
hydrodynamic rotational instability, should show high ampli- 
tude variability as indicated by numerical simulations. Magneto- 
Hydrodynamic simulations of radiation pressure dominated ac- 
cretion disks use however the shearing box geometry, and hence 
they do not provide information about the nature of the viscous 
instability (Hirose et. al. 2009). 

It is interesting to study what effect the viscous stress fluc- 
tuations will have on the overall stability of radiation pressure 
dominated disks. If the viscous fluctuation amplitude is small, a 
linear analysis would naturally show that the unstable mode and 
viscous fluctuations would be uncoupled and the variability will 
simply be their sum. In other words, in the linear regime, the per- 
turbations would grow due to the unstable mode and would not 
be affected by the viscous fluctuations. However, when the vis- 
cous fluctuation amplitudes are large, there may be non-trivial 
and non-linear coupling between the two, perhaps leading to in- 
teresting results. This may especially be true when the systems 
evolves due to the instability and the viscous time-scale on which 
the fluctuations occur changes in the disk. Thus, for a proper un- 
derstanding of these effects, we undertake here a numerical com- 
putation to simulate the time-dependent behavior of a standard 
accretion disk with viscous fluctuations. 

The article is organized as follows. In Section 2. we de- 
scribe our model for the viscosity fluctuations, giving also the 
references to a series of previous work where our radiation pres- 
sure instability numerical model was presented and developed. 
In Section 3., we present the results, showing first (Sec. 3.1.) the 
case of stable accretion disk models with viscosity fluctuations, 
and second (Sec. 3.2) the case of unstable models with fluctua- 
tions. In Section 4, we discuss the results and their implications. 
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2. The Numerical Scheme 

We use the 1(+1)-D hydrodynamical code, developed and de- 
scribed first by Janiuk et al. (2002) for the time-dependent evo- 
lution of an accretion disk. Depending on the mean accretion 
rate, the disk may become unstable due to the radiation pres- 
sure domination. The code has been generalized to incorporate 
two phase flows taking into account mass exchange between the 
disk and an overlying corona (Janiuk & Czerny 2005) and the 
possibility of a fraction of the power being transmitted to the 
a jet/outflow (e.g. Kunert-Bajraszewska et al. 2010) As in the 
standard model, the code assumes that the flow is Keplerian but 
allows for the possibility of energy advection to smaller radii. 

We adopt here the kinematic viscosity of Shakura & Sunyaev 
(1973), as less or equal to the sound speed times the disk vertical 
scaleheight, via the dimensionless proportionality parameter a. 
We incorporate viscous fluctuations in the disk by considering a 
time-varying a parameter given by (see Lyubarskii 1997): 



a(r,t)^aQ[l+/3(r,t)] 



(3) 



where ao is a standard constant value on the order of 0.01 - 
0. 1, and the small factor j8(r, t) accounts for the fluctuating part, 
whose time variation is assumed to follow the Markov chain pro- 
cess: 



/?„ = boUn 

where bo is a parameter and 
u„ - -0.5m„_i + e„. 



(4) 



The factor 0.5 in the above equation ensures that the Markov 
chain does not diverge for large n. In general, the modulus of this 
factor should be less than 1 .0, and its value governs the number 
of series steps between the peaks, A^, as well as the amplitude 
of the oscillations, varu. For factor 0.5, we have A^ « 2.6 and 
varu - Q.lSvare (see e.g. King et al. 2004 and Janiuk & Czerny 
2007). The random variable ?„ is generated with a uniform dis- 
tribution between -1 and 1. The timescale of the fluctuation is 
given by the local viscous time-scale at that radius r expressed 
in terms of the Schwarzschild radius 2GM/c^, as: 



Tvisc = -7^(t7)V('") 
ail H 



(6) 



Here H is the height, Q the Keplerian angular velocity and 
/(r) = 1 - VJ/r is the boundary condition to assure that the 
torque vanishes at the last stable orbit, which also implies that 
the density in the disk goes to zero at its inner edge. The coher- 
ent length over which a varies is assumed to be the local height 
H, which is calculated from the total pressure as H - cJO.. 

The computations are numerically expensive and were per- 
formed with the message passing interface technique, on multi- 
CPU Linux computer cluster machines. We had to restrict the 
number of radius bins to be ~ 200. For typical parameters used 
in this work, this translates to a bin size Ar ~ Hoir)/3, where 
Ho{r) is the height of the steady state disk. This radial grid is 
kept fixed through the time-evolution i.e. we do not adaptively 
change the radial grid. However, in the model, a variability is 
taken to be coherent over a region ~ H, where H is the height at 
a given time. We take this into account by incorporating a mul- 
tiplicative factor on/? of y/Ho/H, which is typically small. This 
factor increases/decreases the viscous variability when the co- 
herent length scale is smaller/larger than the grid size bin. This 
scheme is similar to one used by Mayer & Pringle (2006) to 
model magnetic dynamo evolution. 
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Fig. 1. Lightcurves of models for stable disks where the viscous 
stress scales only with the gas pressure. The green (red) line 
shows the model with viscous fluctuation amplitude bo - 0.6 
(0.99), and the thick black line shows the model with bo = 0.0. 
The mean viscosity parameter ao = 0.1, the mean accretion rate 
is in = 0.33 and the black hole mass is 3 x lO^Mp,. 



3. Results 



(5) 3.1. Stable accretion disk models with viscosity fluctuations 



We first study the effect of viscous fluctuations on stable accre- 
tion disks. Ideally, one should analyze accretion disks where the 
accretion is low, such that the disk is not radiation pressure dom- 
inated. However, for such disks the height H « r and hence the 
numerical computation would require a large number of radial 
bins, to correctly incorporate the viscous fluctuations. Hence, 
we consider disks where the viscous stress scales only with gas 
pressure i.e. Tr^ - aPgas- 

In the presence of viscous fluctuations the luminosity of 
the disk varies, as shown in in Figure [1] where we plot the 
lightcurves of a disk around a 3 x 10** M© black hole with accre- 
tion rate m - 0.33 in Eddington units and q-q = 0.1. As expected, 
as the flickering amplitude bo increases so does the variability of 
the luminosity. It is more illustrative to see the power density 
of the luminosity fluctuations which are shown in Figure |2] For 
small variations in the viscous fluctuations, the variation in the 
luminosity can be estimated to have a power spectrum, roughly 
oc 1 // (Lyubarskii 1997). This is qualitatively consistent with the 
power spectra shown in Figure|2] especially for small bo, i.e. the 
power spectrum is roughly cc 1/f such that /x power is nearly 
independent of frequency. For larger values of bo there is a sig- 
nificant increase of power in low frequency. The reason for this 
behavior is more evident when we plot the light curves for a disk 
around a IOMq black hole in Figure [3] and the corresponding 
power spectra in Figure |4] The lightcurves show that for large 
bo apart from increased stochastic variability, there are also long 
time scale "jumps" in the luminosity. These low frequency vari- 
ations which occur in the viscous time-scale of the outer radii, 
may indeed be due to outer boundary conditions imposed in the 
numerical scheme. The outer disk radius for the computations 
has been taken to be ~ 1000 Schwarzschild radii. Increasing the 
outer radii, incorporates a large cost in computation time espe- 
cially since the total simulation time has to be increased to be 
longer than the viscous time-scale of the outer radius. 
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Fig. 2. Power density spectra of the luminosity fluctuations for 
stable disks where the viscous stress scales only with the gas 
pressure. The models have flickering parameters of bo = 0.6, 
bo - 0.99, bo - 1.5 and bo - 2.0, as shown with green, red, 
magenta, and blue points, respectively. The mean viscosity pa- 
rameter ao - 0.1, the mean accretion rate is m - 0.33 and the 
black hole mass is 3 x IO^Mq. 




Fig. 4. Power density spectra of the luminosity fluctuations for 
stable disks where the viscous stress scales only with the gas 
pressure. The models have flickering parameters of bo - 0.1, 
bo = 0.3, bo = 0.6, bo - 0.99 and bo - 2.0, as shown with cyan, 
magenta, green, red and blue points, respectively. The mean vis- 
cosity parameter ao - 0.1, the mean accretion rate is m - 0.33 
and the black hole mass is IQMq. 
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Fig. 3. Lightcurves of models for stable disks where the viscous 
stress scales only with the gas pressure with flickering parameter 
bo = 0.3, 0.6 and 0.99, marked with magenta, green and red 
lines, respectively, and the thick black line shows the model with 
bo = 0.0 The mean viscosity parameter ao - 0.1, the mean 
accretion rate is m - 0.33 and the black hole mass is IOMq. 



While the low frequency variability at large bo may be an 
artifact of the outer boundary conditions, we note that it will 
not effect the results of our primary study, which is the eff'ect of 
fluctuations on unstable disks. Here the instability time-scale is 
the viscous time-scale of the inner accretion disk region which 
is significantly shorter than the viscous time-scale of the outer 
radius. 

Thus, in general, our numerically computed luminosity vari- 
ations of the viscous fluctuations are consistent with that ex- 
pected from analytical estimates that predict the power spectra 



3.2. Unstable accretion disk models with viscosity 
fluctuations 

Having verified that the results of our numerical code are gener- 
ally consistent with that expected for stable accretion disk where 
analytical results are possible, we consider the more physically 
consistent situation where the viscous stress scales with the total 
pressure. Our primary result is illustrated in Figure For low 
values of viscous fluctuation amplitude bo - 0.66, the disk is 
unstable and shows periodic very large amplitude outbursts as 
expected. 

It can be seen from the plot that the strong fluctuations in 
the luminsoity occurs when the flux level of the source is high. 
This is because during such "high" flux states, the disk is thicker 
and hence the stochastic viscous fluctuations occur over a larger 
coherent range. Moreover, the viscous time-scale is shorter and 
hence the effects of viscous fluctuations are more pronounced. 

Increasing bo - 0.99 causes the oscillations to become ir- 
regular and the duration of the bursts also changes. Further in- 
crease of bo - 2.0 decreases the amplitude of the bursts and the 
lightcurve becomes more "chaotic" in nature. Finally for a large 
bo - 2.5, the lightcurve after a transient beginning phase, be- 
comes stochastic and there is no evidence for any cyclic bursts. 

The time averaged accretion rate in all the simulations shown 
in Figure|5]is equal to the external accretion rate supplied to the 
disk {in - 0.33). When the radiation pressure instability oper- 
ates, the disk luminosity cannot be constant, corresponding to 
this average acretion rate, but switches between the low and high 
luminosities, corresponding to the stages when the accretion rate 
onto the black hole is either much higher or much lower than the 
average value. When the instability is suppressed by the viscous 
fluctuations, there are no longer any 'low' or 'high' states. 

For stellar mass black holes, we find similar results with 
the disk showing oscillatory behavior for low values of bo and 
stochastic behavior for higher values (see Figure |6). The disks 
stabilize at a lower value of /jq ~ 1 as compared to bo > 1 .5 
for AGN disks. This maybe related to the fact that at the same 
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Eddington rate, an AGN disk is more radiation pressure domi- 
nated than that for a stellar mass black hole. 

We also consider the effect of power being transferred to a 
jet. Specifically, we parametrize the jet (cf. Janiuk & Czerny 
20 11) by including an additional term in the energy equation 



Q 



cP, 



rad 



kI. 



(1 - ?/jet) 



(7) 



where the fraction ?7jet is assumed to depend on on time and dis- 
tance from the black hole via the local accretion rate: 



'7jet = 1 - 



1 



1 + AjetOT(r, t)^ ' 



(8) 



with Ajet being a parameter of the model. The physical meaning 
of this parametrization is that at the Eddington accretion limit 
there would be equipartition between the jet and disk radiation 
for Ajet = 1- We adopt a value of Ajet = 5 or = 10 , which means 
a moderately strong jet. Without viscous fluctuations, even for 
such moderately strong jets the disk undergoes cyclic oscilla- 
tions, but with a reduced amplitude. The inclusion of compara- 
tively small viscous fluctuations bo - 0.3 stabilizes the disk for 
both values of A^^,,. Thus a moderately strong jet decreases the 
amount of viscous fluctuation required for stability. 

In Table [T] we present the summary of the models we have 
calculated. In summary, large amplitude viscous fluctuations 
stabilize the disk against the radiation pressure instability. The 
strength of the viscous fluctuation required to stabilize the disk 
is slightly lower for stellar mass black hole systems than AGNs 
and decreases when there is an effective jet in the system. The 
Table also shows the normalized rm.s variability of the lumi- 
nosity. Note that for locally stable disks (i.e. when t,^ oc Pgas), 
the r.m.s increases with b,,, while for locally unstable disks(i.e. 
when Tr^ oc Ptot), it generally decreases. This is because the vis- 
cous fluctuations stabilize the disk and high amplitude oscilla- 
tions disappear. 

We defer a detailed comparison of these results with the ob- 
servations of a particular source for future work. This maybe 
done for the various states of the microquasar GRS 1915-1-105 
(Belloni et al. 2000). Another good candidate maybe the recently 
studied source IGR J 1709 1-3624, which also exhibits the so- 
called 'heartbeat' states for luminosity variations in some spec- 
tral states, as well as stochastic fluctuations (Altamirano et al. 
201 1). The detailed modeling would require a proper determina- 
tion of the Eddington ratio as well as an estimation of the fraction 
of energy going into jet power from radio observations. 



4. Summary 

Our numerical simulations show that large amplitude viscous 
fluctuations can stabilize radiation pressure dominated standard 
accretion disks. This may explain why even at luminosities close 
to the Eddington, black hole systems exhibit spectral signatures 
of a standard accretion disk which do not show large amplitude 
cyclic behavior Instead, they show stochastic variation which 
maybe due to viscous fluctuations which stabilize the disk. 

At low luminosities, the disk would be gas pressure domi- 
nated and hence intrinsically stable. Here, viscous fluctuations 
may produce variability in the luminosity. However, since it is 
likely that the coherent length of the viscous variability scales as 
the height of the disk which is much smaller than the radius for 
low accretion rate systems, the net effect on the global luminos- 
ity should be small. In other words, since H/R « 1 for these 
sources, the viscous fluctuations at different radii will incoher- 
ently add to give a diminishing effect on the total luminosity of 
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Fig. 5. Lightcurves in the 5 models of accretion disk oscilla- 
tions due to the radiation pressure instability limit-cycle. The 
results are shown for the alpha flickering models with parame- 
ters bo = 0.0, bo = 0.6, bo = 0.99, bo = 2 and bo = 2.5, from 
top to bottom, respectively (the quiescent period between thew 
outbursts in bo - 0.0 model is about 0.6 Myrs). The model pa- 
rameters are: black hole mass of 3 x 10** Mq, accretion rate of 
m - 0.33 in Eddington units, ao = 0.1. The model includes no 
jet. 



the disk. This may explain why the standard disk component in 
some low luminosity systems are found to be steady with little 
variability. 

While, in this work, we consider some general characteris- 
tics of Galactic black holes and AGNs, a more detailed com- 
parison with data of one or two specific sources may provide 
both qualitative and quantitative insights. In particular, the time- 
varying spectra and not just the integrated luminosity from these 
numerical simulations should be compared with data from some 
chosen sources. The model should be confronted with both the 
variability as a function of energy as well as frequency depen- 
dent time-lags. 

In the light of the results of this work, it will also be inter- 
esting to re-examine the models for GRS 19 15 H- 105 for cases 
when the system shows oscillations as well as when it does 
not. Detailed comparison between the expected and observed 
lightcurves, may reveal whether the viscous variability ampli- 
tude changes between the different states or it is the fraction of 
energy going into the jet. The study may finally reveal the distin- 
guishing feature of the source as compared to other black hole 
systems. 

Finally, a more comprehensive numerical analysis needs to 
be undertaken with different boundary conditions and assump- 
tions, that may identify the crucial conditions necessary for the 
unstable disk to be stabilized by viscous fluctuations. An adap- 
tive radial grid will be in particular useful. The rotating disk may 
also give rise to non-trivial effects, which can be studied only by 
a full 3-D hydrodynamic code. Indeed, what seems to be impor- 
tant is that the non-linear effect of large amplitude perturbations. 
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Table 1. Summary of the time-dependent disk models. The ac- 
cretion rate in Eddington units is m = 0.33 and mean viscosity is 
ffo = 0.1 for each case. rm.s. is the normaUzed root mean square 
variability of the luminosity. 
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Fig. 6. Lightcurves in the 5 models of accretion disk oscillations 
due to the radiation pressure instability limit-cycle. The results 
are shown for the alpha flickering models with parameters /jq - 
0.0, /7() = 0.3, bo = 0.6, bo = 0.99 and bo = 2.0, from top 
to bottom, respectively. The model parameters are: black hole 
mass of IOM0, accretion rate of m = 0.33 in Eddington units, 
ao = 0.1. The model includes no jet. 

which can be captured only by numerical simulations, can give 
rise to non-trivial results such as stabilization of a system which 
is locally highly unstable. 
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